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(N ' 
O 

q ' The measurement of the trilinear couplings A in the part of the Lagrangian which breaks 

. supersymmetry softly will be a difficult experimental task. In this report the heavy Higgs 

decays H,A^> f\fi to stau pairs are investigated for measuring the stau trilinear coupling 
A T . Based on detailed simulations of signal and backgrounds for a specific reference point 
in future high luminosity e + e~ linear collider experiments, it is concluded that the param- 
Q-i! eter A T can be determined with a precision at the 10% level in the region of moderate to 

J|h \ large tan (5. 



1 Introduction 
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The couplings between fermionic matter fields and Higgs fields differ from those of the scalar 
matter fields once supersymmetry is broken, see e.g. Ref.[l]. In theories based on soft super- 
symmetry breaking the scalar-Higgs Yukawa couplings are modified multiplicatively by the A 
parameters which, in parallel to the fermion-Higgs Yukawa couplings, are inter-generational 
matrices. In accordance with bounds on flavor-changing couplings the A parameters are gener- 
ally assumed to be diagonal and three parameters, A t , At and A T , are introduced for the third 
generation. 

By definition, the A parameters come with the Yukawa couplings which are of the size of 
the fermion masses. Therefore they cannot be measured in general directly except for the third 
generation. Since they couple Higgs fields with scalar L-fields and i?-fields, they become effec- 
tive in two ways: (i) They contribute to the off-diagonal elements in the scalar mass matrices, 
and to the mixing of L- and -R-states; and (ii) They give rise to mixed scalar L and R decay 
final states of the heavy scalar and pseudoscalar Higgs bosons. 



In the scalar stop sector the off-diagonal mass matrix element is given by m t (A t — \i cot j3). 
For moderate to large tan j3 the second term is suppressed and A t can be determined quite 
accurately by measuring the mixing effects in the stop mass spectrum and the stop mixing angle 
in e + e~ annihilation to stop pairs [2]. In heavy Higgs decays, on the other hand, A t is shielded 
by the potentially much larger term //tan/3 and Higgs decays to stop pairs, if kinematically 
allowed at all, are less suited for measuring the stop trilinear parameter. 

The situation is reversed in the down sector, i.e. for staus. While in the stau system A T 
is shielded by the term //tan/3 in the mass matrix [3], the A T parameter is enhanced by the 
coefficient tan (3 in the couplings of the heavy scalar and pseudoscalar Higgs bosons to mixed 
pairs of tl and ?r fields. Heavy Higgs decays are therefore promising channels for measuring 
the stau trilinear parameter A T . 

The expressions for the partial decay widths become especially transparent in the limit 
where (i) the heavy Higgs boson masses are large [decoupling limit], (ii) tan/? is large, and 
(iii) the LR mixing is small. In this limit the decay widths of the scalar and pseudoscalar Higgs 
bosons to mixed pairs fif 2 = t+t^ + ff f 2 + are given by 

r(H,A^nf 2 ) * gg^i/sOMan/g + Ai) 2 

4a/27t m HiA 

where A accounts for the pase space suppression in the usual form. The couplings of the scalar 
Higgs boson H to diagonal pairs of L- and i?-fields are suppressed by coefficients m T /A T and 
m z / (A T tan P) which both are small in the limit we are considering. The coupling of the 
pseudoscalar Higgs boson A to diagonal pairs vanishes in CP-invariant theories. 

Using the partial widths for Higgs decays to tau pairs, 

T(H,A^tt) ~ ( fI^L miiA tan 2 /?, (2) 

4V27T 

the decay widths to stau pairs may be normalized by the decays to tau pairs: 

T(H, A ^fif 2 ) _ Al/2 (Ar + ficotp) 2 (3) 



F(H, A^tt) m 



H,A 



If the normalization is chosen alternatively by the dominating bb final states, the ratio of the 
widths is reduced by a coefficient m T /3ml: 



T(H, A — » fif 2 ) _ 1/2 m 2 (At + \i cot (3) 
T(H, A —> bb) ~ 3m 2 m H 



xl/2 m r l^r + /i CQipj 

A — 2 — 9. • ^' 



In any case, for moderate to large tan (3 and A T of the same order as /t, the size of the branching 
ratio of the heavy Higgs bosons to mixed LR stau pairs is essentially set by A 2 T . Thus, for 
sufficiently large A T the measurement of these branching ratios provides a valuable instrument 
for measuring A T . 



2 



2 Properties of the Higgs system 



The qualitative arguments presented above appear strong enough to perform a quantitative anal- 
ysis in order to prove this method to be useful for measuring A T in practice. For this purpose 
we adopt the mSUGRA reference point SPSla' defined for the SPA Project [4]. It is closely 
related to the standard reference point SPSla, yet with a cold dark matter density in accordance 
with the WMAP measurement. 

The mSUGRA parameters are defined as M = 70 GeV, M 1/2 = 250 GeV, A = 
—300 GeV, tan/3 = 10 and sign [i = +. Extrapolation to the electroweak scale generates 
the Lagrangian parameters A T = —445 GeV and fx = 403 GeV, thus \A T \ ^> //cot/3 holds 
indeed. The masses and branching ratios of the supersymmetric particles relevant to the present 
analysis are summarized in table 1 . 



Particle 


Mass [GeV] 


Decay B 


Decay B 


H° 


431.1 


t-t+ 0.075 
bb 0.683 
it 0.053 
rff^ 0.014 
ffff 0.031 
f 2 ~f 2 + 0.003 


XiX\ o.on 
X°iX°2 0.040 
X2X2 0.023 
XtXi 0.056 


A 


431.0 


t-t + 0.055 
bb 0.505 
ti 0.103 

f^ff 0.035 


X°iX°i 0.011 
0.055 
X^X^ 0.063 
XtXi 0.170 


x\ 


97.8 






X\ 


184.4 


Tfr* 0.564 
eje* 0.024 
0.026 


i/ r z/ r 0.155 
v e v e 0.115 
0.115 


xi 


184.2 


t+v t 0.519 


z> T r+ 0.189 
v e e + 0.138 
i> M /i + 0.138 


fi 


107.4 


XiT~ 1.000 






195.3 


XiT~ 0.869 
X° 2 r- 0.046 


XiV T 0.086 


V T 


170.7 


X°u T 1.000 





Table 1: Masses and branching ratios of heavy Higgs bosons, light gauginos and third gener- 
ation sleptons in the SPSla' scenario [4]. The Higgs decays are calculated with the program 
FeynHiggs 2.2.10 [5] 

At a linear collider with an energy y^s of about 1 TeV heavy Higgs boson production 
e + e _ — > HA, see [6], will clearly be kinematically accessible for this reference point [7]. 
The measurement of their decay modes, however, will confront several problems: 

- Due to their mass degeneracy the decays of H and A cannot be resolved experimentally. 
Thus one can only determine the branching ratios for the sum of both Higgs bosons. 
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- The energy spectra of the final r decay products reflect only weakly the energy of the 
primary particles, which is gradually softened during cascade decays involving massive 
invisible particles like neutralinos or sneutrinos. It is therefore extremely difficult to dis- 
criminate f\ from f 2 decays. Instead, only the sum of all fjfj decay modes will be deter- 
mined. 

- As a consequence of the (moderately) large value tan/3 = 10 the neutralino an d 
chargino xf decays lead preferentially to final states involving r leptons. Abundant multi- 
tau signatures constitute a severe background to all channels involving SUSY particles, 
in particular to the decays of interest H,A^ fifj. 

The strategy to determine the H, A decay modes and branching ratios is the tagging of one 
Higgs particle by its decay into a pair of bb jets and the analysis of the recoiling system: 

e+e" -> HA -> bbX . (5) 

The decay modes and event topologies under investigation are 

= f\T\+T\Tz + f%Tz — > T + T~$, (6) 

X^o = xlxl + xlxl - r + r-$, (7) 

X^- = xtxl - T+T~$, (8) 
and the reference decay modes are 

X TT = T + T~ , (9) 

X b i = bb^ jet jet. (10) 

The particle content of the three supersymmetric final states (6) - (8) is identical. In order 
to distinguish these channels it will be assumed that the masses of the primary and of all the 
secondary SUSY particles are known well enough so that the resulting decay topologies and 
r ± spectra can be reliably modeled and simulated. This knowledge is important in order to 
determine the branching ratios of the various decay modes from their relative contributions to 
the 'observable' data distributions. This assumption is quite natural as the measurement of the 
A parameters is certainly a second-generation task. Details of the event generation are presented 
in the appendix A. 

The cross sections for HA pair production [6] assuming common scalar and pseudoscalar 
Higgs masses are shown in figure 1 for the three center of mass energies ^fs = 0.8 TeV, 
1.0 TeV and 1.2 TeV. The remaining parameters are taken from the reference point SPSla' 
for illustration. The present study is representative and based on 10,000 HA events which, 
for scenario SPSla', may be accumulated with a cross section of 1.8 fb at 1 TeV or 3.9 fb at 
1.2 TeV, respectively. The results may be easily scaled to lower statistics event samples without 
losing their significance. 



3 Experimental analysis 

In this section the analyses of reaction (5) with the Higgs decay modes (6) - (10) will be de- 
scribed in detail. As mentioned above the channels involving supersymmetric particles, Xf^, 
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Figure 1: Cross section for e + e~ — > HA production as a function of the common H, A mass 
at y/s = 0.8 TeV, 1.0 TeV and 1.2 TeV. The curves include e ± beam polarizations of V e - = 
±0.9 and V e + = ±0.6, as well as QED radiation and beamstrahlung effects 

X^o^o and X-+~-, lead to the same final state and the topologies do not allow their separa- 
tion on an event-by-event basis. Rather a statistical analysis will be applied to determine their 
branching fractions. The decays into Standard Model particles, X TT and X b i, can be efficiently 
isolated and will be used for normalization. The results will be given in terms of combined 
branching ratios, defined as B bix = B(H -> 66) B(A -> X) + B(A -> 66) B(H -> X). 

It should be noted that the results for the background Higgs decays to charginos and neutrali- 
nos can probably be predicted at the time of the stau analyses. They depend only on parameters 
which can be measured in the chargino/neutralino sector itself at earlier times. This way the 
experimental results of the Higgs decays to charginos and neutralinos can be compared with 
theoretical predictions. 

3.1 Signal channel e + e~ — »■ HA — ► bb r + T~Ifi 

The topology is characteristic for all Higgs decays into supersymmetric particles. The criteria 
listed in table 2 are chosen in order to optimize the acceptance for HA — > bb f»fj — > bb t + t~ Ifi 
decays. 

The criteria (1) - (3) provide a very efficient selection of HA — > bbX events by tagging one 
Higgs particle via its resonant decay into a pair of b quark jets [see discussion in appendix A 
and figure 6]. The good energy resolution allows the reliable transformation into the rest frame 
of the recoil system X which is identified as the second Higgs particle. 

The criteria (4) - (6) select SUSY decays into secondary r's plus large missing energy. 
The last cut (7) removes direct decays into rr pairs, which are back-to-back in the Higgs rest 
frame. The properties of the various decay modes are displayed in the left panels of figure 2, 
where normalized distributions of the visible tau energy and di-tau mass are shown. It is a 
common feature of both spectra that the dominant contributions come from XiXi peaks at low 
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Selection criteria 


Constraint 


1 two identified b jets 




2 b jet energy 


100 GeV <E b < 400 GeV 


3 bb invariant mass 


™>h,a — 30 GeV < m bb < m H ^ + 30 GeV 


recoil mass against bb 


m H,A — 30 GeV < m recoi i < m Hj A + 90 GeV 


4 two oppositely charged r candidates 




5 visible r energy 


2.5 GeV <E T < 200 GeV 


visible tt energy 


E TT < 250 GeV 


6 missing energy 


250 GeV < E < 550 GeV 


7 acollinearity angle in Higgs rest frame 


C T > 10° 



Table 2: Event selection criteria for the signal reaction HA — > bb fjfj — > bb r + r Ifi 
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Figure 2: Spectra from e + e~ — * iJA — > bb t + t~ $ decays of a) visible tau energy E* in the 
H, A rest frame; b) di-tau mass m TT . Left: normalized distributions; right: fitted contributions 
of individual channels fjfj, xl,X2 an d XiXi t0 me observable signal. mSUGRA scenario SPSla' 
at Vs = 1 TeV 
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values, while the spectra from fjfj and XiX.2 extend towards higher values. The separation 
between the decay modes improves when correlations between both particles are exploited, e.g. 
the invariant mass m TT . Also notice that the shapes of the distributions from X1X2 an ^ X2X2 
are barely distinguishable, thus only the sum of all neutralino decays, labeled x?X2' w iU be 
investigated. 

The overall H, A — > fifj efficiency is ~ 43%. However, there are still large contributions 
from Higgs decays into charginos (~ 37%) and neutralinos (~ 23%), both of which have higher 
combined branching ratios, see table 1 . 

The distributions from the complete simulation of the visible r energy in the Higgs rest 
frame E* and the di-tau mass m TT are shown in the right panels of figure 2. The contributions 
from the individual decay modes, X f . f ., X^o and X-+-- of eqs. (6) - (8), are summed up and 
fitted to reproduce the data of figure 2. The analyses of the observables E* and m TT emphasize 
different characteristics but lead to consistent results. In all fits (including observables not 
shown) the chargino contribution can be determined in a stable manner whereas the stau and 
neutralino parts are strongly correlated. 

The fit results are displayed in the spectra of figure 2. The relative rates, acceptances and 
the extracted combined branching ratios B bbx are summarized in table 3. 



e+e" -> bbX 


fbhX 6 bbX 


&bbX 


HA -> bb Tifj 
bb XiX°2 
bbxtXi 


0.186 ±0.041 0.428 
0.292 ±0.052 0.228 
0.516 ±0.036 0.372 


0.049 ±0.011 
0.135 ±0.024 
0.146 ±0.010 


HA -> bbrr 


0.515 


0.075 ±0.004 


HA -> bbbb 


0.630 


0.345 ±0.007 



Table 3: Expected accuracies on the determination of Higgs decays HA — > bbX. Listed are the 
analyzed event samples, the fitted contributions of decay modes fjx x , the detection efficiencies 
e b i x , and the combined branching ratios B b i x - The results are based on 10,000 HA decays in 
the SPS la' scenario 

3.2 Reference channels e + e~ —*■ HA — > bb tt and HA — > bb bb 

e + e~ — > HA — > bb t + t~ The selection of HA — ► bb t + t~ decays is complementary 
to the analysis of the previous SUSY decays. The basic criteria (1) - (4) of table 2 for bb 
and tt identification are applied. However, the r energy spectra are harder and both t's are 
emitted back-to-back in the Higgs rest frame, leading to the following cuts: (5) visible r energy 
5 < E T < 400 GeV, tt energy E TT < 500 GeV; (6) no missing energy requirement; (7) 
acollinearity angle in the Higgs rest frame £* r < 10°. 

The reconstructed spectrum of the visible r energy E* in the Higgs rest frame, shown in fig- 
ure 3, is fairly flat and extends up to the energy of the primary, undecayed r lepton. The overall 
detection efficiency is high, see table 3. A combined branching ratio of B b i TT = 0.075 ± 0.004 
can be obtained, where only statistical uncertainties are given. The analysis may be further 
improved by an overconstrained kinematical fit. Exploiting energy-momentum conservation 
and approximating the r directions by the directions of the decay products and treating the r 
energies as free parameters, allows one to construct 2 constraints (2-C fit), see [8]. 
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Figure 3: Spectrum of the visible r energy E* in the H, A rest frame from e + e" 
bb t + t~ decays. mSUGRA scenario SPSla' aty/s = l TeV 



HA 



e+e — > HA — > bb bb The selection of HA — > bb bb events is straightforward by applying 
the same criteria (1) - (3) of table 2 to another pair of bb jets. The four jets are then combined 



(i) 

bb 



m 



(2) 
bb 



such as to construct two bb systems with invariant masses closest to each other, m 
Again, the selection efficiency is high. The energy distribution of the b jets in the Higgs rest 
frame, displayed in figure 4, exhibits a clear signal of a narrow peak at half the Higgs mass. The 
combined branching ratio for the decay mode X bb can be determined with a statistical accuracy 
of B b - bbi = 0.345 ± 0.007, for details see table 3. 




Figure 4: Spectrum of 6-jet energy E b in the H, A rest frame from e + e" 
decays. mSUGRA scenario SPSla' at -s/s = 1 TeV 



HA 



bb bb 



The measurement only provides information on the product of cross section times branching 
ratio. In order to extract the decay rates the HA production cross section has to be calculated 
accurately, which in turn requires a precise knowledge of the Higgs masses. The b jet energy 
distribution of figure 4 (or equivalently the bb mass spectrum similar to figure 6a) can be used 
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to determine the H, A masses with an accuracy of 5m H ^ ~ 0.15 GeV. This error can be 
reduced further by applying kinematic fitting techniques [8]. In fact such a procedure allows 
the selection of a very clean HA event sample with low background, so that relaxing of the 6 
quark identification criteria may be envisaged. 



4 Interpretation and conclusions 



The expected results for the combined branching ratios of Higgs decay modes (6) - (10) are 
summarized in table 3. The double ratios are experimentally determined as B bbfif . / B bbbb = 
0.142 ± 0.032 and B bbf . f ./B bbrT = 0.653 ± 0.147. Their relations to the partial decay widths 
can be written as 1 



B h 



B 



bbbb 



B(H -> 66) B(A -> fif 2 ) + B{A -> bb) B(H -> 



B(H - 



bb) B(A 



r(i4 -> 66) r(# -> 66) 



66) 



(11) 



B 



bb : 



B(H -> 66) £(A -> fif 2 ) + £(A -> 66) B(# -> f^) 



66 1 



66) £(A 
^ flf 2 ) 



rr 



hS(A 



► 66) 



rr 



T(A -> rr) (1 



1 

2 



r(i4 



rir 2 J 



+ 



T(H - 



rr) (1 + 1/r) 



rr 



_ T(A^rr) r(/f 
In the second double ratio the two terms in the denominator have been identified, i.e. 

T(A -> 66) r(ff -> rr) 



F(H -> 66) T(A -> rr) 



(12) 



(13) 



which is expected to hold with high accuracy. 

These double ratios are proportional to (A T + fx cot P) 2 ~ in the decoupling limit for 
large A T and large tan j3 when Li? decays dominate over the diagonal LL and RR decays and 
mixing can be neglected. For the parameters chosen in this study, however, we must include 
corrections from LR mixing of the particles and the diagonal LL and RR decays. 

The pseudoscalar Higgs boson A couples to off-diagonal fif 2 pairs with the same amplitude 
as to tlTr pairs so that no mixing corrections need be applied. In contrast, the coupling of the 
scalar Higgs boson H to off-diagonal stau pairs is modified by the mixing parameter cos 29 f 
and, moreover, H decays include also contributions from the genuine diagonal couplings tlTl 
and trTr of the order m T /A T and m z /{A T tan (3) with respect to the leading off-diagonal am- 
plitudes. Since the mixing parameter 



sin 29 f 



2 



(A T — fi tan/3) 



(14) 



'As mentioned earlier, the pseudoscalar Higgs boson A decays only to off-diagonal T1T2 pairs in CP-invariant 
theories while the scalar Higgs boson H can decay to all combinations of stau pairs 
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Figure 5: Double ratios of the combined branching ratio B HA ^ bbf . fj normalized to B HA ^ bbbb 
(left) and B HA ^ bbTT (right) as a function of A T . The lower curves show the contributions from 
diagonal pairs f 2 f 2 (blue) and fifi (magenta). The horizontal (green) lines indicate the expected 
experimental accuracy based on 10,000 HA decays in scenario SPSla' 

involves A T itself, the dependence of the decay amplitudes on A T is modified and not linear 
anymore. As a result, the binomial character of the partial decay widths in A T is changed 
asymmetrically. 

The dependence of the double ratios of the decay widths, c/eqs. (11), (12) and (3), on A T 
including the subleading LR mixing effects are calculated using FeynHiggs [5] and are dis- 
played in figure 5 together with the expected experimental accuracies. There are two possible 
solutions for A T at —450 GeV and +350 GeV, which, however, can be distinguished experi- 
mentally because they correspond to different f mixing configurations. The mixing parameter 
sin 29 f differs by ~ 20 % for the two solutions. Since the mixing can be obtained with an accu- 
racy of a few percent from measurements of the fi mass and the fifi production cross section 
at e + e~ colliders (see [9, 3] for scenarios with similar parameters), this additional information 
is sufficient to single out the negative A T solution. 

From the simulation of heavy Higgs decays into supersymmetric and SM particles one ob- 
tains for the trilinear coupling 

A T = -450 ±50 GeV 

as an ab-initio determination of this soft SUSY breaking parameter for an event sample of 
10,000 HA Higgs pairs. 

Uncertainties in the theoretical predictions and parameters, eqs (3) - (4) and their analogues 
for diagonal decays, are expected to be negligible at the level of achievable experimental ac- 
curacies. Theoretical calculations are under control at the per-cent level when all the one-loop 
corrections in the f/r and Higgs sectors are included [10]. It is interesting to note that the 
parameter tan/5 can be controlled internally within the same analysis of HA — > bbbb de- 
cays. The tan (3 dependence of the combined branching ratios can be expressed as B bbbb = 
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1/[(1 + c A /tan 2 /5)(l + c^/tan 2 /?)] with coefficients c A ^ 100 and c H ~ 50 for A and if 
decays, respectively, in the reference point considered. From the measurement quoted in ta- 
ble 3 one expects a precision of S tan (3 ~ 0.15. The parameter /x can be measured in chargino 
production within a few per-cent. Both uncertainties result in a shift of the trilinear coupling 
of at most 5A T < 1 GeV, far below the anticipated experimental error. These estimates are 
confirmed by a combined analysis of SUSY parameters based on measurements of many SUSY 
production processes at the ILC and LHC [1 1]. 

The direct determination of the trilinear coupling analyzed in the present report may be com- 
pared with other methods which make use of higher order corrections affected by the parameter 
A T . A global analysis by means of Fittino [11] provides a combination of X T = A T — [i tan (3 = 
-4450 ± 30 GeV, together with tan (3 = 10.0 ± 0.1 and \i = 400.4 ± 1.3 GeV. However, ap- 
plication of this indirect method is possible a priori only in scenarios in which the degrees of 
freedom are specified in toto when the virtual loop corrections are included and if all theoretical 
uncertainties are under proper control. In contrast, the proposal described in the present paper 
is a robust leading order analysis. 

5 Summary 

While the trilinear stop-Higgs coupling A t can be measured fairly easily by evaluating the stop 
masses and the mixing angle, this task is much more demanding for the trilinear coupling A T 
in the stau sector since these couplings come with the masses of the quarks and leptons. Nev- 
ertheless, we have demonstrated in this report that the measurement of A T is possible in scalar 
and pseudoscalar Higgs boson H, A decays. Large luminosities at the e + e~ linear collider 
ILC would be required, however, to achieve an accuracy of about 10%. Though the measure- 
ment is difficult, this direct determination based on tree-level processes is necessary before the 
determination through indirect effects based on quantum corrections can be trusted with high 
confidence. 

After the stop trilinear coupling A t will be determined, the measurement of at least one 
additional trilinear parameter is required to investigate universality properties of these parame- 
ters, for instance, as implemented in minimal supergravity. The A T measurements are therefore 
important ingredients for reconstructing the underlying physics scenario [12]. 
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FeynHiggs program. The work of SYC was supported in part by the Korea Research Foun- 
dation Grant (KRF-2002-041-C00081) and in part by KOSEF through CHEP at Kyungpook 
National University. PMZ thanks the German Science Foundation DFG and Stanford Univer- 
sity for partial financial support during an extended visit of SLAC where part of this work was 
carried out. 

A Event generation 

Events are generated with the program PYTHIA 6.3 [13] which includes initial and final state 
QED radiation as well as beamstrahlung a la ClRCE [14]. The decays of r leptons are treated by 
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TAUOLA [15]. The detector simulation is based on the detector proposed in the Tesla TDR [7] 
and implemented in the Monte Carlo program Simdet 4.02 [16]. The main detector features 
are excellent particle identification and measurement of charged and neutral particles for a polar 
angle acceptance 9 (n — 9) > 125 mrad. 

In the analysis the reconstructed b jets and r candidates are required to be within the accep- 
tance of | cos 9\ < 0.95, while Higgs boson are produced centrally ~ sin 2 9. The identification 
of HA — ► bbX events is provided by the good jet energy flow measurement with a resolution 
of a/E = 0.3/ y/E(GeV). This is illustrated in figure 6, where the di-jet mass and the mass 
recoiling against the bb system are shown. Both distributions are fairly narrow and peak at the 
Higgs masses. The recoil mass spectrum is slightly wider and extends towards large values due 
to radiative effects. 

For the identification and reconstruction of r candidates, a narrow jet with invariant mass 
m T < 2.5 GeV is required which contains one charged particle plus possibly additional photons 
or three charged particles. In general the leptonic 3-body decays r — > ev e v T (17.8%), r — > 
^ v ii v T (17.4%) are less sensitive to the primary r energy than the hadronic decays r — > ixv T 
(11.1%), r — > 7r ± 7r°z/ T (25.4%) and r -> 'k ± -k^-k~v t + -k ± -k°'k g v t (19.4%). All decay modes 
are used in the analysis, except of ee and jiji pairs. 
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Figure 6: Spectra from e + e — > HA — > bb X decays of a) bb di-jet mass m^; b) recoil mass 
rrirecoii against bb system. mSUGRA scenario SPSla' at y/s = 1 TeV 

Since the decay rates of interest are of the order of a few percent, to be further degraded 
by efficiency losses, the case study is based on a high statistics sample assuming a production 
rate of Nha = &HA ■ C = 10, 000 events. The total cross section of HA production amounts 
to a ha — 1-8 fb at — 1 TeV, including e ± beam polarization, QED radiation and beam- 
strahlung, see figure 1 . The results of the present study may be easily transferred to any other 
energy or reference point once the parameters are specified. 

The characteristic event signatures HA — > bbrr E, i.e. two energetic b jets forming a high 
mass resonant state plus two r leptons plus (possibly large) missing energy, are very clean. Any 
background from QCD processes qq(g), WW or ZZ production is estimated to be small and is 
therefore neglected. 
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